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o\ ; 

. Some exact results are derived relating quantities computable by the so-called damage spreading 

' method and thermodynamic functions for the N-color Ashkin-TeUer model. The results are valid 

for any ergodic dynamics. Since we restrict our analysis to the ferromagnetic case the results are 
also valid for any translational invariant lattice. The derived relations should be used in order to 
determine numerically the N-color Ashkin-TeUer critical exponents with better accuracy and less 
computational efforts than standard Monte Carlo simulations. 



Since its introduction, the damage spreading (DS) method has become a powerful tool in the study of phase 
transitions and critical phenomena [ll, Q • Basically it consists in following the time evolution of two initially identical 
^ copies of a given system, both subject to the same dynamical rules and to the same thermal noise, after the introduction 
^ of a small perturbation (called damage) in one of them at a given time. This method presents some advantages on the 
^ traditional Monte Carlo (MC) method where the time evolution of a single copy of the system is investigated For 
I instance, in the DS method the fluctuations are substantially reduced as compared to the conventional MC method, 
"j^ \ A recent review on the DS method can be found in ^] . 

■ In statistical mechanics, the most important question when performing a simulation is how to obtain the equilibrium 
^ ; ' thermal properties of the system under investigation. Thus it is very important to establish exact relations relating 
1 measurable quantities obtained in a DS simulation to the thermodynamical properties of the system. Such relations 
Ch ' were first obtained for the ferromagnetic Ising model in a square lattice and the numerical results showed a prominent 
I ■ reduction of fluctuations and finite-size effects In the wake of those findings the exact relations between damage 
' O ■ spreading computable quantities such as magnetization and spin-spin correlations functions and thermal equilibrium 
C \ quantities has been obtained for several systems like, for instance, the Potts [6*1, Ashkin- Teller discrete N- vector 
Q ' Hi (Na^Np) [^], and spin-1 Q models. These relations are valid for any translationally invariant system which 
evolves in time under any ergodic dynamics. Recently, such relations were applied to numerically investigate the 
Potts Tol, two-color Ashin- Teller [ll| and a spin-1 [9] models. In all these cases the computed critical exponents were 
(N in aggreement to the exact known results to within four decimal places, revealing the accuracy of the DS method. 

The Ashkin-TeUer (AT) model was introduced to investigate cooperative phenomena in quaternary alloys . It 
was shown to be equivalent to a two-layer Ising model with a four-spin coupling between the two layers |13j . In 
ly-^ , two dimension the AT model can be mapped onto a staggered eight-vertex model at the critical point and shows 



> 
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■ non-universal critical behavior along a self-dual line where the exponents vary continously In general, despite 



its simplicy the two-color AT model displays a wide variety of critical and multicritical phenomena already in two 
dimension, as has been shown in studies based on du ality arguments [l^ , real-s pac e renormalization group [3 [3 , 



(yQ mean- field renormalization-group approachs [19|, |20|, |2l| , finite-size scaling [22|, |23|, |2J| , and conventional Monte Carlo 
O • simulations [2^, [II . 

IL^ ' The N-color Ashkin- Teller (N-AT) model consists of N Ising models coupled pairwise through a four-spin interaction 
[29} . When = 2 it becomes the usual AT model discussed above (also known as two-color AT model). It was argued 
that for > 2 the model has a first-order transition as long as the four-spin coupling is ferromagnetic [2^ . In [29j | the 



H ' phase diagram was investigated by mean-field analysis and Monte Carlo simulation for N = 3 and by several other 
- - I tecniques for A^ > 2 [s^, [lU, HI, [H, [H, [H, [H, [13] ■ Recently there is a renewed interest in this model mainly due 



to the richness of its phase diagram and its relation to other systems [38|, |39|, |40|, |41|. However, it still lacks a lot of 
work in order to fully understand the main features of the N-AT for A^ > 2 in any dimension. We believe that the 
DS simulations are quite suitable to give us more information about this model. 

The purpose of the present note is to determine exact relations envolving thermodynamic quantities and some 
specific combinations of damages for the N-color Ashkin-TeUer model. The model is defined by the Hamiltonian 
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where — ±1, a and & labels different Ising spins (distinguished by their "color") at a site i and {ij) indicates that 
the sum is performed over all distinct pairs of sites on a given lattice. For A'^ = 1 we recover the Ising model, whereas 
for iV = 2 we have the usual Ashkin- Teller model. In what follows, our treatment is valid for any value of N , and for 
( J2, J4) in the ferromagnetic range in order to take into account the translational invariance of the lattice. 
According to the usual mean-field procedure we can write a single-site effective Hamiltonian 



N N 

meff = E + 2 E ^'^>'^ (2) 

where /i° = K2^j^^{(Jj) and h'^^ — K^'^-_^^{a'^(j^). In the above equation (• ■ •) denotes thermal average and where 
we have introduced the new variables Kn — Jn/ksT [n = 2,4). The structure of the effective fields acting on af and 
(jf<j\ implies that the order parameters are given by 



m'^^K), a=l,.--,Ar, (3) 

and 

M'^^^{ato\), a^b, (4) 
while the correlation functions are expressed by 



ij 



(5) 



and 



r^" = - «a,^) (a^a,^) , a ^ b. (6) 

To implement the damage spreading technique it is convenient to introduce the following binary variables 



- ^(i + O, (7) 

nf = lil + ala',), a^b. (8) 

Let us consider two configurations (referred to as A and B) on a regular, translational invariant, lattice evolving 
in time according to the same (and ergodic) dynamics such as Metropolis, heat-bath or Glauber [3|. We define four 
different types of damage combinations: 



(1) n^(A) = l and nnB)=0, (9) 

(2) n^(A)=0 and n^(i3) = 1, (10) 

(3) nf(^) = l and nf(B) = 0, a^b, (11) 

(4) nf(^) = and nf(B) = l, a^b. (12) 

It is known that after a long time the system eventually reach the thermal equilibrium. In this regime, the above- 
defined damages occur with probabilities given, respectively, by: 



p, = {[l-Ut{A)]unB))t, 
P3 = {Uf{A)[l-Uf{B)]), 
P4 = {[l-Uf{B)]ufiB))t 



(13) 
(14) 
(15) 
(16) 
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where (• ■ ■)t means time average over the trajectory followed by the copies of the system in their phase space. In what 
follows it is convenient to introduce differences between such probabilities: 

F^p,-p2^{Uf{A))t-{Il^iB))t , a=l,---,7V, (17) 

and 

G = p3 - P4 = (nf (A))t - (nf (B))i , a^b. (18) 

The next step in our analysis consists in imposing some constraints in the temporal evolution. In the present case 
it turns out that there are four distinct such possibilities: 

(ei) copy A evolves without any constraint, while copy B is restricted to the constraint that for an arbitrarily fixed 
site, say i 0, ng(B) = 0, for a 1, ■ • • , TV. 

(62) copy A evolves with ng(A) — 1 and copy B with ng(i?) = 0, for a = 1, ■ ■ ■ , N. 

(63) copy A evolves without any constraint, while copy B is restricted to with Il^^{B) =0, a^h. 

(64) copy A is subjected to ng''(A) — 1, while copy B is restricted to ng''(i?) =0, a^h. 
Ergodicity implies, for the evolution (ei), that 

(n?(A)), = (nn (19) 

and, with the help of conditional probability, 

(n?(s))* = . (20) 

Eqns (|19l20p relates, in a definitive way, the dynamical and thermal averages represented, respectively, by (• • ■)t 
and (••■). Substitution of those equations into p?|) gives us 

. ,21, 

which, with the help of ([3]), ([5]) and ([7]), can be expressed as 

^(^1) ^ on^°' ■ (22) 
2(1 — m°) 



For the evolution [e-i) we have, in a similar way. 



and 



from which follows 



m{A))t = (23) 



(11" fi - n")) 



4 



In order to write the final Eqs. (|22p and ([25)1 we have used the fact that m° = as a consequence of the lattice 
translation invariance. The results expressed by these equations means that at first we compute numerically F[ei) and 
F{^2) (quantities related to DS). Then, we may determine the magnetization and the two-point correlation functions 
(thermal equilibrium quantities) as 

= 2^ - 1, (26) 



and 



F(e2) 



rS.^4|;M(F(e2)-nei)), (27) 
^ (62) 



for a = 1,2, •••,7V. 

For the evolution (ea) we find 



and 



{ILf{A))t - (nf ) (28) 



(nf (i?))t = ^ \ ^ ^' . (29) 



1 - (ng^) ' 



for a ^ b. 

Thus, from (g]), (O and (O we get 



Finally, evolution (54) implies that 



G(e3) = ^^l, a ^ 6. (30) 
^ ' 2(1 - M"'') ^ ^ 



(nf (A))t = ^ " ^ (31) 



and 



(nf(i-nf)) 

for a ^ b 



{nf{B))t = , 7' , (32) 



Therefore, in terms of M'^'' and Fg^ we have 

G(e4) = YT^|^, (33) 
Thus, having determined G^e^) and G(e4) in a numerical simulation, we can use (|30p and (|33[) to compute 

j^ab ^ _ 1 (34) 

G(e4) 



and 



rS.'' = 4^(G(e4)-G(e3)). (35) 
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FIG. 1: The correlation function r'i'(r) versus r. The slowest decay was found for T/Tc = 1.002, in which case the plot of 
logj^Q[r^^'(r)] versus logjQ r is presented in the inset, leading to the estimate 77 = 02521 ± OflOSO. 



Since F{ei), F{e2), Gi^s) and G{e4) are computed from selected combinations of damages, the relations expressed 
by ([2H-I2S1) and ([501 - which are exact for any translationally invariant lattice at all temperatures, allow us to 
determine numerically the thermal properties of the TV-color Ashkin- Teller model. 

As an example we considered the ferromagnetic 3-color Ashkin- Teller model on a square lattice of linear size i = 60 
and with J4 = ODIJ2. Periodic boundary condition were applied and we computed the two-site correlation functions 



r'^'Hr) = -Y.r'.f (a = 1,2,3), (36) 

i(r) 

measured with respect to the central site located at (L/2, L/2). The two-site correlation functions are computed 
from (j27p . As the equilibrium was approached we found that F(i)(r) = F(2)(r) = F(3)(r) as expected from the 
permutation symmetry between the spin variables presented by the Hamiltonian. Initially a copy denoted by A with 
all spin variables trl"^ (A) = 1 is left to evolve according to a MC run during a time teq — I x 10'* which is sufficient 
to the system get close enough to thermal equilibrium. Then a second copy denoted by B is created by replicating 
copy A and introducing small modifications, the initial damage, on it. Both copy are left to evolve under the same 
ergodic dynamical rules during a time teq = 18 x 10^ and subjected to specific boundary conditions corresponding 
to the possibilities (ei) and (62) above mentioned. In addition to that, each simulation were performed for M = 20 
different samples in order to reduce inherent statistical fluctuations for a given temperature. To locate the transition 
point we search for the temperature ratios T/Tc at which the functions F^")(r) show the slowest decay as a function 
of r. At criticality we expect that, for large r. 



rW(r)^r-'' (a = 1,2,3). (37) 

In Fig. 1 we present the behavior of F^^^ (r) for several values of T/Tc (here Tc is the usual Ising critical temperature 
at the decoupling point J4 = 0). It is noted that the slowest decay occurs for T/Tc = lfl02. The inset in Fig. 1 
shows a log-log plot of F(*)(r) versus r, at the observed critical point, from which we obtain rj = 02521 ± 0D050. 
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Therefore, for the 3-color AT model we found an evidence of a continuous transition of Ising type for J4/ J2 = ODl. The 
complete phase diagram for this particular case can be determined by the present method, as long as the transitions 
are continuous. We hope to address to this problem in a future work. 

Recently published results [§, [l^, [ll| in other nontrivial situations have established the accuracy of the present 
approach to obtain critical properties of classical spin models in translationally invariant systems. We thus have 
presented news exact relations conecting some damage-spreading functions to thermal equilibrium properties for the 
ferromagnetic N-color Ashkin- Teller model and we hope that such relations will become a useful guide to high-precision 
numerical investigations. 
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